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Abstract 

A simple lower bound to the capacity of the discrete-time Poisson channel with average energy 
e s is derived. The rate | log(l + e s ) is shown to be the generalized mutual information of a modified 
minimum-distance decoder, when the input follows a gamma distribution of parameter 1/2 and mean 

I. Introduction 

Consider a memoryless discrete-time whose output Y is distributed according to a Poisson distri- 
bution of parameter X, the channel input. By construction, the output is a non-negative integer, and 
the input a non-negative real number. The channel transition probability W(y\x) is thus given by 

W{y\x) = e~ x — (1) 

This model, the discrete-time Poisson (DTP) channel, appears often in the analysis of optical com- 
munication channels. In this case, one can identify the input with a signal energy and the output with 
an integer number of quanta of energy. 

Let Px{x) denote the probability density function of the channel input. We assume that the input 
energy is constrained, i. e. E[X] < e s , where E[-] denotes the expectation operator and e s is the average 
energy. Random variables are denoted by capital letters, and their realizations by small letters. 

An exact formula for the capacity C(e s ) of the DTP channel is not known. Recently, Lapidoth and 
Moser |T], derived the following lower bound 

C(„)> log ((l + I) 1+ V)-( 1 + ^). (2) 

Observe that this bound diverges for vanishing e s . Capacity is given in nats and the logarithms are 
in base e. 
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A closed-form expression for the mutual information I(X; Y) achieved by an input with a gamma 
distribution of parameter v was derived by Martinez in 0, namely 

/<*n=/i«.-d w ^ d " 




(y + e s (l — u)Y ) 1 — u J logu 

+ (e a + i/)log^±^+e a (#/ + l)-l), (3) 

where ip(y) is Euler's digamma function. For v = 1/2, numerical evaluation of the mutual information 
gives a rate which would seem to exceed \ log(l + e s ) for all values of e s . In this paper, we prove that 
the rate | log(l + e s ) is indeed achievable by this input distribution. The analysis uses a suboptimum 
minimum-distance decoder, similar in spirit to Lapidoth's analysis of nearest neighbor decoding 0. 

II. Main Result 

Let the input X follow a gamma distribution of parameter 1/2 and mean e s , that is, 

Pjr(x) = -=L=e-£. (4) 

V27T£ S 2 

This choice led to good lower and upper bounds in [2 and respectively. 

We consider a maximum-metric decoder; the codeword metric is given by the product of symbol 
metrics q(x, y) over all channel uses. The optimum maximum-likelihood decoder, for which q(x, y) = 
W(y\x), is somewhat unwieldy to analyze (Eq. (0 gives the exact mutual information). We consider 
instead a symbol decoding metric of the form 

q (x,y) = e- ax -^, (5) 

where a = 1 + The reasons for this choice of a will be apparent later. 

Clearly, the decoder is unchanged if we replace the symbol metric q(x, y) by a symbol distance 
d(x, y) = — log q(x, y), and select the codeword with smallest total distance, summed over all channel 
uses. This alternative formulation is reminiscent of minimum-distance, or nearest-neighbor decoding. 
Indeed, the metric in Eq. ® is equivalent to a minimum-distance decoder which uses the distance 

d( x , y)= ^ V ~ ^ X)2 = V - + ax- 2y^. (6) 

x x 

The term —ly^fa is common to all symbols x and can be removed, since it does not affect the 
decision. 

For a = l, the distance in Eq. © naturally arises from a Gaussian approximation to the channel 
output, whereby the channel output is modeled as a Gaussian random variable of mean x and variance 
x. This approximation is suggested by the fact that a Poisson random variable of mean x approaches 
a Gaussian random variable of mean and variance x for large x. 
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Minimum-distance decoders were considered by Lapidoth in his analysis of additive non- 
Gaussian-noise channels. For our channel model, even though noise is neither additive (it is signal- 
dependent), nor Gaussian, similar techniques to the ones used in Q can be applied. More specifically, 
since we have a mismatched decoder, we determine the generalized mutual information [4]. For a 
given decoding metric q(x, y) and a positive number s, it can be proved [4] that the following rate 
— the generalized mutual information — is achievable 

q(X,Y) s 



Xgmi(s) = E 



log' 



(7) 



/ P x {x)q(x',Y) s dx' 

The expectation is carried out according to Px(x)W(y\x). This quantity is obviously a lower bound 
to the channel capacity. 
Our main result is 

Theorem 1. In the discrete-time Poisson channel with average signal energy e s , the rate \ log(l + e s ) 
is achievable. 

This rate is reminiscent of the capacity of a real-value Gaussian channel with average signal-to-noise 
ratio e s . Similarly to the situation in this channel, the rate is achieved by a form of minimum-distance 
decoding. Differently, the input follows a gamma distribution, rather than a Gaussian. 

Proof: We evaluate the generalized mutual information Igmi(s) for an input distributed according 
to the gamma density, in Eq. (01). First, we evaluate the expectation in the denominator |5, Eq. 3.471- 
15] 







e 2es x e 

— dx' = - (8) 

V2vre s x' y 1 + 2ae s s 

1}Q 



Further, using the expression of the first two moments of the Poisson distribution, namel) 

Y / W{y\x)y = x, ^W{y\x)y 2 = x 2 + x, (9) 

y y 

together with the input constraint J Px{x)xdx = e s , we can explicitly carry out the expectation in 

'The moment generating function of a Poisson random variable of mean x is readily computed to be e x ^ e The first 
two moments are the first two derivatives, evaluated at t = 0. 
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Eq. ©, 



-J P x {x)Y^W{y\x)\og(^J P X (x')q(x' \y) s dx^j dx (10) 
i j P x {x) W(y\x) (^-ax - dx 

-J P x (x)J2w(y\x) (-yJ 2sil + * a£sS) -\ogy/l + 2ae s s\ dx (11) 



= -s ((a + l)e, + 1) + y/2e 8 s(l + 2ae s s) + - log(l + 2ae s s). (12) 

Choosing s = 7 — Jf* , rrrrr i the first two summands cancel out. And for a = 1 + — we 
have that 2as = 1, and therefore 

/gmi(s) = -log(l + e s ). (13) 



The same rate, | log(l + e s ), is also achievable by a decoder with a = 1. In this case, we have to 
replace the generalized mutual information by the alternative expression I L m El, given by 

a(X)q(X,Y) s 



'LM 



E 



log- 



(14) 



/ P x (x)a(x')q(x',Y) s dx' _ 
As for Igmi, s is a non-negative number; a(x) is a weighting function. Setting a(x) = e~~ x we have 
that I L m is given by Eq. (fTTT) . thus proving the achievability. 

The bound provided in this paper is simpler and tighter than Eq. It would be interesting to 
extend Theorem[T]to channel models Y = S(X) + Z, where S(X) corresponds to the case considered 
here and Z is some additive noise Z, with a Poisson or a geometric distribution. A different input 
distribution and another modified decoding metric are likely required for either case. 
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